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The phase diagram and critical behavior of the spin-1 Baxter- Wu model with a crystal field in 
two dimensions is explored by renormalization group, conventional finite-size scaling and conformal 
invariance techniques. We found that the phase diagram of this model is qualitatively the same 
as that of the dilute 4-states Potts model, presenting a multi-critical point for a finite value of the 
crystal field, in disagreement with previous work based on finite-size calculations. However, our 
results indicate that the critical exponents vary continuously along the second-order transition line, 
differently from the expected behavior of the dilute 4-states Potts model. 
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I. INTRODUCTION 

The Ising modeli was the first non-trivial model ex- 
actly solvable in two dimensions which exhibits sponta- 
neous symmetry breaking. It became the most popular 
ferromagnetic model in statistical mechanics and even 
today is the object of several studies in other contexts 
like random systems^ The dynamics of the Ising model 
is described by the Hamiltonian 

Hi = -J ^ s i s 3' (!) 

<i,j> 

where the sum is over all nearest neighbors and the clas- 
sical spin variables Sj = ±1 are attached at each site i of 
the lattice. In middle 60's, Blume and Capel? proposed 
an extension of the Hamiltonian (TIJ to study first-order 
magnetic phase transitions. Their Hamiltonian is given 

by 

h bc = - j s ^ + A Y, s *> ( 2 ) 

<i,j> i 

where A plays the role of a crystal field and in this case 
the variables Si are classical spin-1 variables taking the 
values Si = — 1, 0, 1. 

It is well established that for dimensions d > 2 the 
Blume-Capel model (Eq. <|2j|) presents a phase diagram 
with ordered ferromagnetic and disordered paramagnetic 
phases separated by a transition line which changes from 
a second-order character (Ising type) to a first-order one 
at a tricritical point (see Ref. |j and references there in) . 
More specifically, in two dimensions, the machinery com- 
ing from conformal invariance^ indicates that at this 
multi-critical point the long-range fluctuations are gov- 
erned by a conformal field theory with central charge 
c = 7/104i& In this case, all the critical exponents and 
the whole operator content of the model were obtained 3 
(see also Ref. 0). The generalization to higher spin S 
of this model has also been studiedi 4 ^' 1 ^ 11 ' 12 In particu- 
lar, results of mean field theoryjii conformal invariancep 



and Monte Carlo simulations^ predict different phase 
diagrams for integer or half-odd-integer spins, in contra- 
diction with results based on real space renormalization 
groups Aifi Recently, the universality at a double critical 
endpoint in the two-dimensional spin-3/2 Blume-Capel 
model has been analyzed and it was shown that it be- 
longs indeed to the same universality class as the critical 
line»A£ 

Another simple model exactly solvable in two dimen- 
sions exhibiting spontaneous symmetry breaking is the 
Baxter- Wu modelm^^^ This model is defined on a tri- 
angular lattice by the three-spin interaction Hamiltonian 

Hbw = -J X! SiS i Sk ' ( 3 ) 

<ijk> 

where the sum extends over all elementary triangles of 
the lattice and Si = ±1 are Ising variables located at 
the sites. This model is self-dualA£ with the same criti- 
cal temperature as that of the Ising model on a square 
lattice. The critical behavior of the Baxter- Wu model is 
governed by a conformal field theory with central charge 
c = and its leading exponent o 14 ^^ a = 2/3, 

v = 2/3 and rj — 1/4 are the same as those of the 4- 
states Potts models 

In analogy with the Blume-Capel model, in this paper 
we consider the Baxter- Wu model in the presence of a 
crystal field. The Hamiltonian of the spin-1 Baxter- Wu 
model with a crystal field is given by 

H = -J J2 s * s r s k + A Y, S ^ ( 4 ) 

<ijk> i 

where the classical spin variables Sj, defined in a trian- 
gular lattice, take the values s, = —1,0,1. Note that 
when A — > — oo only the configurations with Sj = ±1 are 
allowed, and we recover the pure Baxter- Wu model. 

Since we have in the Baxter- Wu model with a crystal 
field the same kind of competition between the ordered 
((s) ^ 0) and the disordered phases ((s) = 0) (which 



2 



is mediated by the crystal field) as in the Blume-Capel 
model, we may expect for both models a similar phase 
diagram, but with different critical behavior. This kind 
of competition also appears in the dilute g-states Potts 
model, 1 - It is well known that the Baxter-Wu model and 
the 4-states Potts model have the same critical exponents 
(see for example Ref. 14). Since the dilution in the 4- 
states Potts model has the same effect as the crystal field 
in the Baxter-Wu model, we may expect the critical be- 
havior of both models to be the same. Some previous 
calculations of the phase diagram of both models have 
been reported in the literature. Nienhuis et al.f£ based 
on a renormalization-group study, indicate that for the 
dilute 4-states Potts model the phase diagram is similar 
as that of the Blume-Capel model, i.e., there is a tran- 
sition line which changes from a second-order character 
to a first-order one at a multi-critical point. In this case, 
however, the critical behavior is governed by only one 
fixed point, giving along the second-order line the same 
exponents as that of the pure 4-states Potts model. On 
the other hand, Kinzel et al.^ using finite-size meth- 
ods, conjectured a different kind of phase diagram for 
the Baxter-Wu model in the presence of a crystal field. 
These authors interpreted the changes of the estimated 
thermal exponents yt along the transition line as a sig- 
nal that a second-order transition should happen only for 
A —>■ — oo (the pure Baxter-Wu model). 

A careful study of the pure Baxter-Wu model has also 
been done by exploring its Bethe-ansatz solutionpiii^*i^ 
Using the consequences of conformal invariance, it has 
been show n 1 ^ 1 ^ that in the absence of dilution (A — > 
— oo ) not only the leading critical exponents of the 
Baxter-Wu model and the 4-states Potts model are iden- 
tical, but the whole operator content of the models co- 
incides as well. Moreover, the masses of the field the- 
ory describing the thermal and magnetic perturbations 
of both models are also identical^ Since the reported 
effects of the dilution in both systems are different, we 
decided in this paper to study the effect of dilution (or 
crystal field) in the Baxter-Wu model in two dimensions 
using the machinery of conformal invariance and renor- 
malization group techniques. 

This paper is organized as follows. In the next section 
we present the phase diagram obtained through the mean 
field renormalization group approach. In section III we 
present the transfer matrix of the model, the relations 
used in our finite-size studies, as well as the results for 
the phase diagram. We then close the paper in section 
IV with a summary and conclusions. 



II. MEAN FIELD RENORMALIZATION GROUP 

The mean field renormalization group (MFRG)2ii2£ is 
a powerful phenomenological approach which can pro- 
vide quite good results in the general study of critical 
phenomena.^ It is based on the comparison of the or- 
der parameter for different finite lattices in the presence 



of symmetry breaking boundary conditions. For a finite 
cluster of N spins, and considering the parameters of the 
Hamiltonian @J, one first computes the magnetization 
per spin m,ff(K, D, b), where K = f3J, D = /3A and b is 
the boundary field, with /3 — l/fcgT and fcs the Boltz- 
mann constant. As the boundary field is assumed to be 
very small one has 

m N (K,D,b) = f N (K,D)b. (5) 

In its simplest version, the MFRG considers two dif- 
ferent clusters of N' < N spins and assumes that the 
magnetizations scale as 

m N >(K', D', b') = £ d - yH m N {K, D, b), (6) 

where £ — [N/N') 1 ' d is the scaling factor, d is the dimen- 
sion of the lattice and yn the magnetic exponent. With 
the same relation for the boundary fields, namely 

b'=l d - yH b, (7) 

and taking into account expansion Q one has 

f N (K,D) = f N ,(K',D'), (8) 

which is independent of any scaling factor and is viewed 
as a renormalization recursion relation from which one 
gets fixed point solutions K' = K = K c and estimates 
of correlation length critical exponents v in the subset 
D 1 = D. 

Another version of the approach considers three differ- 
ent clusters of N" , N' and N spins (in increasing order) 
together with the scaling © in such a way that one gets 

f N ,(K',D')b' =t d 1 - yH f N (K 1 D)b 1 (9) 



f N „(K",D")b" = £ d - yH f N ,(K\D')b\ (10) 

where the scaling factors are l\ — (N/N 1 ) 1 ^ and £2 = 
(N'/N") 1 / d . Imposing now that the boundary fields scale 
not as bulk magnetizations above, but as the surface field 
we obtain 

b' = l\ h °b, b"=t y 2 h °b', (11) 

where yi ls is the surface critical exponent. Eqs. 10-dJ 
are now the renormalization recursion relation between 
the interaction parameters of the system. The exponent 
d — yh — yhs is then determined self-consistently by im- 
posing further that Eqs. © and 11 Oil possess the same 
fixed point K" = K' = K = K c for an invariant subset 
D" = D' = D (for further details and a comparison be- 
tween these two approaches see Refs. and . This 
version of the method is referred as surface bulk MFRG 
(SBMFRG). 

Before applying the MFRG to the spin-1 Hamiltonian, 
defined by Eq. Q, it is worthwhile to measure the effi- 
ciency of the method by first treating the pure Baxter-Wu 
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model (Eq. where exact results are available. Since 
the model exhibits three different ferrimagnetic phases at 
low temperatures^! it is supposed, at first sight, that the 
size of the finite blocks must be such that they will suit- 
ably accommodate them in an equivalent way through- 
out the lattice. These would dramatically restrict the 
size of the clusters. However, we note that by comput- 
ing the sub-lattice magnetizations itia, tub and mc, by 
taking three different boundary fields i>A, &b and be, we 
obtain the same equations as by considering a homoge- 
neous cluster where tua = tub = nic and bA = &s = be- 
This is not surprising, bearing in mind that the ferromag- 
netic phase is also coexisting with the other three ferri- 
magnetic ones at low temperatures. This means that, 
for practical purposes, we can consider only the ferro- 
magnetic arrangement, allowing us to take blocks of any 
number of spins. This leads a huge simplification in the 
numerical acquisition of the functions Jn{K). For this 
particular system, the best clusters are those which pre- 
serve the symmetry of the original lattice. They are made 
of symmetric triangles of N = 6, 10, 15, 21, 28 spins. Al- 
though for N = 6 we can obtain analytical expressions for 
me(K, b) and fe(K), for N > 10 all the quantities have to 
be computed numerically. In Table|U we present the crit- 
ical temperature and critical exponent obtained accord- 
ing to the usual MFRG and also from the SBMFRG. It is 
also possible through the MFRG get an extrapolation^ 
Note that, for the present case, the extrapolated value is 
not so close to the exact result. This is, in fact, expected 

Table I: Results for the pure Baxter-Wu model according to 
the MFRG and SBMFRG approaches. 



N'-N 


K c 


V 




MFRG 




6-10 


0.28667 


2.78203 


10-15 


0.36565 


1.89869 


15-21 


0.39335 


1.55875 


21-28 


0.41077 


1.37178 


extrapolated 


0.43768 


0.69995 


exact 


0.44069 


2/3 




SBMFRG 




6-10-15 


0.58138 


1.06875 


10-15-21 


0.48090 


1.08061 


15-21-28 


0.47832 


0.99573 



since the MFRG does not reproduce the exact value as 
the size of the lattices go to infinity^ Although the SBM- 
FRG doeS)22 the present clusters are still too small for 
us to get a reasonable extrapolation. However, as can be 
seen from Table [IJ the temperatures are getting closer to 
the exact value as the size of the clusters increase. Note 
that, while the MFRG approaches the expected result 
from below, the SBMFRG does it from above. A com- 
mon feature of the present renormalization group is still a 



worse estimate of the critical exponent when compared to 
the critical temperature. Nevertheless, in general, a rea- 
sonable picture of the criticality of the system is achieved 
from the approach. 

We now proceed to the study of the spin-1 model with 
crystal field anisotropy. The same symmetry arguments 
also apply for this model. However, due to computer 
time, we were here limited to block sizes N = 6, 10, 15. 
The phase diagram in the 5 = A/J versus kT/J plane is 
depicted in Fig. ^according to both procedures, as well 
as from the finite-size scaling (FSS) procedure of the next 




Figure 1: Phase diagram in the S versus kT/J plane for the 
spin-1 Baxter-Wu model in the presence of a crystal field. 
Continuum lines are second-order phase transitions and the 
dashed line is a first-order transition (see next section). The 
dot represent the multi-critical point. The results are ac- 
cording to the MFRG with two cluster, SBMFRG with three 
clusters, FSS (see next section) and MF approximation with 
one cluster. The sizes of the clusters used are also presented. 

section, and the usual mean field (MF) approximation. 
The latter approach, is obtained by assuming mjv = b 
in Eq. ©, resulting in Jn = 1. Except for the FSS, all 
the lines terminate at some point which is identified as 
the multi-critical point. The first-order transition lines 
are not possible to be obtained neither from the MFRG 
nor from the SBMFRG. Comparing with the FSS result, 
discussed in the next section, we can see that the MFRG 
overestimates the critical temperature while the surface 
bulk version underestimates it. This is just what hap- 
pens in the pure case, as shown in Table Q] when com- 
pared to the exact critical value. ( nao entendi o q vc 
quer dizer... In all cases, the critical temperature of the 
spin- 1/2 model is obtained in the limit 5 — > — oo. pode 
tirar esta frase? Ja' nao foi dita na introducao!) The 
estimate of the multi-critical point is given in Table [H] 
The SBMFRG and MFRG give critical exponents that 
vary along the second-order critical line. However, this 
fact, may be just an artifact of these aproaches, since 
there is only one renormalization recursion relation^ In 
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the next section, we analize in the context of confomal 
invariance the possibility of the critical exponents vary 
along the second-order line. 



Table II: Position of the multi-critical point for the curves 
shown in Fig. U 



N' — N 




i>t 




Mb Klj 




6-10 


1.1816 


2.1523 


10-15 


0.9330 


1.8462 




SBMFRG 




6-10-15 


0.6408 


1.5835 




MF 




6 


0.3133 


0.8141 


10 


0.3539 


1.0328 


15 


0.4513 


1.1902 




FSS 




3-6-9 


1.2225 


1.3089 



III. 



FINITE-SIZE SCALING AND CONFORMAL 
INVARIANCE 



The row-to-row transfer matrix T of the Hamil- 
tonian (| in a triangular lattice, with horizontal 
width N, has dimension 3^ x 3 , Its coefficients 
(s[, s' N \T\si, sn) are the Boltzmann weights gen- 
erated by the spin configurations {si, s^} and 
{s[, s' N } of adjacent rows. If we consider periodic 
boundary condition in the horizontal direction, the trans- 
fer matrix can be written as 



(Si, S N \f\s[, 



■J s n) 



N 



[]exp [t-hi+rs'tist + s' i+1 ) - St^s 2 ] , (12) 

with t = ksT/J {ks is the Boltzmann constant) and 
5 = D/J. 

The finite-size behavior of the eigenvalues of T 
(A (iV) > Ai(iV), ...) can be used to determine the criti- 
cal line and the critical exponentsASi^ The critical line 
(t c (S)) is evaluated by extrapolating to the bulk limit 
(N —>■ oo) the sequences i c (<5, N) obtained by solving 



G N {t c )N = G N+3 (t c ){N + 3), N = 3,6,. 



(13) 



where GN{t c ) is the mass gap of H = — lnT and is given 
by 



G N {t c ) = In 



Ao(AQ 
A x (iV) 



The multi-critical points are obtained using a heuris- 
tic method, which has already been proved to be 
effective4^ In this case we have to solve simultaneously 
111 .'it for three different lattice sizes 



G N (t c )N = G N+3 (t c )(N + 3) = 



G N+6 (t c ){N + 6), N = 3,6,. 



(14) 



In Eqs. Ill -it and 11411 we restricted the possible finite 
strip widths to multiples of 3 in order to preserve the 
invariance of the Hamiltonian <@J under the reversal of 
all spins on any two sub-lattices. 

As usual, we expect the model to be conformally invari- 
ant in the region of continuous phase transition. This in- 
variance allows us to infer the critical properties from the 
finite-size corrections of the eigenspectrum at t c £& The 
conformal anomaly c, which labels the universality class 
of critical behavior, can be calculated from the large-A^ 
behavior of the ground-state energy of H = — In T 6 



mA (AQ 
TV 



6N 2 



■o{N~ 



(15) 



where eoo is the the ground-state energy per site in the 
bulk limit and v s — V3/2 is the sound velocity. The scal- 
ing dimensions of operators governing the critical fluctu- 
ations (related to the critical exponents) are evaluated 
from the finite-A^ corrections of the excited states. For 
each primary operator, with dimension x^, in the oper- 
ator algebra of the system, there exists an infinite tower 
of eigenstates of H = — lnT whose energy hi(A^ m ,) and 
momentum Pf n , are given by^ 



lnA t, m 'W lnA o (A0 2irv s 



N N N 2 

P m,m> = -^-{S4> + m-m 



( x 4> 



m - 



(16) 



where to, to' = 0, 1, ... . 

A finite-size estimate for the first-order transition line 
can be obtained by the same procedure done as in a 
recent worki^ At the first-order line of the Baxter- Wu 
model with a crystal field, we have the coexistence of 
five phases, one ordered ferromagnetically, three ordered 
ferrimagnetically and a disordered one. Consequently, 
for a given lattice size we calculate the points where 
the gap corresponding to the fifth eigenvalue has a mini- 
mum. The extrapolation A^ — > oo of these points give us 
an estimate for the first-order transition line. 

In the numerical diagonalization of 11211 we used 
the Lanczos method for non-Hermitian matrices^ We 
also considered the translational symmetry to block- 
diagonalize the transfer matrix. 

In Fig. Q]we show the second-order transition line (con- 
tinuum line), obtained by solving Eq. l(T3jl for lattice 
sizes N — 6. As we can see in this figure the second- 
order transition line also occur for finite values of the 
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crystal field (S ^ — oo), differently of the previous results 
of Kinzel et aL^H where it was conjectured the appear- 
ance of the second-order transition only at 8 — ► — oo. 
For the pure Baxter- Wu model {5 — > — oo) we obtained 
t- 1 = 0.4408842 for N = 6 in Eq. Q3J, which differs only 

0.04% of the exact value t' 1 = ln(v f +1) = 0.440686... . 
For this reason, it is a very good approximation to con- 
sider t c (S) = t 6 c ' 9 (S). 

For the sake of clarity, we present in Table lUll the finite- 
size sequences obtained by solving Eq. l(T3)) for lattice 
sizes N = 3 and N = 6. As we can see from this table 
the convergence of t c are better for 5 > 0, so we expect 
that the estimates of {&) are better than the corre- 
sponding ones in the region 8 < 0. The fast convergence 
with TV indicates that the corrections to finite-size are 
probably given by a power law, like in the pure Baxter- 
Wu model. In this last case the corrections are controlled 
by an operator with dimension w = 4iA& 

Table III: Finite-size data t^' N+3 given by 11311 for the critical 
temperature i c for some values of S. 



N 


-10 


-1 


1 


1.25 


1.3 


3 


2.246498 


1.843818 


1.358399 


1.251121 


1.226940 


6 


2.256769 


1.849705 


1.360144 


1.251529 


1.227005 



model with a crystal field are shown in Table|EI The ex- 
trapolated values of c°° can be obtained from the TV-large 
behavior of c N ' N+3 (5), given by 

c N ' N+3 (S) = c°° -a (N- w + (N + 3)- w ) x 

(N- 2 + (N + 3)- 2 y\ (18) 

where a is a constant and w the dominant dimension 
associated to the operator governing the finite-size cor- 
rections. We determine c°° assuming Eq. l|18Jl exact for 
N = 3 and N — 6 with w given by the pure Baxter- 
Wu model, i.e. w = 4. IS. We see from Table &\ that 
the conformal anomaly is c = 1 along the critical line, 
and apparently even at the multi-critical point. This 
scenario is quite different from that of the Blume-Capel 
model, where the conformal anomaly changes abruptly 
from c = 1/2 at the critical line to c = 7/10 at the tri- 
critical point (see Ref. 5), but its qualitatively similar 
as that of the dilute 4-states Potts model^ where along 
of the critical line the critical exponents as well the con- 
formal anomaly do not change. Since we do not have 
a precise estimate for the multi-critical point, we have 
also calculated c N ' N+3 for several values of t c (6) between 
1.2 < t c (S) < 1.4. We have not seen any abrupt change of 
the conformal anomaly, like in the Blume-Capel model.' 



We have also solved Eq. ll"14t for N ~ 3, which gives 
us an estimate for the multi-critical point. In this case 
we do not have points to extrapolate, but we believe the 
estimate of the multi-critical point, the last point in the 
continuum line in Fig. is not far from the extrapolated 
one. Our estimate for this point is t t — 1.2225 and St = 
1.3089. 

We determined the first order line minimizing the gap 
related with the fifth eigenvalue, as discussed before. The 
dashed line shown in Fig. Q]was obtained in this proce- 
dure considering N = 9, As we can see in this figure, 
the first-order transition line finishes at the multi-critical 
point. 

In the critical regions of the phase transition line (con- 
tinuum curve) the conformal anomaly and the scaling 
dimensions can be calculated exploring the conformal in- 
variance relations II 1511 and llfit . From Eq. Ill oil a possible 
way to extract c is by extrapolating the sequence 



n N,N+3 



12 



lnA 1>0 (iV + 3) lnA li0 (7V) 



7V + 3 



N 



1 



(/V + 3) 2 



1 

N 2 



(17) 



calculated at t c (8). In Eq. ifTTIl A„^(iV) means the nth 
eigenvalue of 111 211 with size N in the sector with momen- 
tum p. Examples of such sequences for the Baxter- Wu 



Table IV: Finite-size estimates c N ' N+3 , given by 117(1 . for the 
conformal anomaly is shown for some values of 5. The ex- 
trapolated results obtained by 118|l is also shown. 



N 


-10 


-1 


1 


1.25 


1.3089 


3 


0.938546 


0.941169 


0.954322 


0.955674 


0.955047 


6 


0.986393 


0.986829 


0.986656 


0.979694 


0.975435 


oo 


1.006 


1.005 


0.999 


0.989 


0.984 



From Eq. Ill (ill the scaling dimensions x(n,p) related 
to the nth (n = 1,2,...) energy in the sector with mo- 
mentum p can be obtained by extrapolating the sequence 



N 



Ai,o(JV) 
ttV3 "* V A «,p(A0 



:ln 



(19) 



In Tablelvlwe show the dimensions x N (l, 0) and x N (2, 0) 
for the Baxter- Wu model with a crystal field. Note that 
x (1, 0) = x 0) due the Eq. E3 an d the fact we choose 
i c (<5) = t®' 9 (5). In the tricritical point the three en- 
tries 2:^(1,0) are same due Eq. El For 5 — » — oo (the 
pure Baxter- Wu model) the scaling dimensions x 9 (l,0) 
and a; 9 (2,0) differ only 1% from the leading dimensions 
x 9 (l,0) = | and a; 9 (2,0) = \ of the pure Baxter- Wu 
model. As the estimate of the critical line is better in 
the region S > and the eigenvalues converge faster with 
the size of the lattice in this region, we believe that our 
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Table V: Finite-size scaling dimensions 2:^(1,0) and x (2,0) 
given by l|19jl for some values of 5. 







N=9 


N=6 


N=3 




5=-W 


0.1236 


0.1236 


0.1195 




5=-l 


0.1223 


0.1223 


0.1184 


x N (l,Q) 


5=1 


0.1113 


0.1113 


0.1093 




5=1.25 


0.1048 


0.1048 


0.1043 




5=1.3089 


0.1026 


0.1026 


0.1026 




5=-W 


0.5059 


0.5147 


0.6057 




5=-l 


0.4897 


0.4976 


0.5756 


x N (2,0) 


5=1 


0.3644 


0.3762 


0.4163 




5=1.25 


0.3023 


0.3207 


0.3613 




5=1.3089 


0.2829 


0.3037 


0.3456 



estimates for the scaling dimensions are better in the re- 
gion 5 > than the corresponding ones for 5 < 0, i.e., the 
estimates x 9 (l, 0) and x 9 (2, 0) must differ by less than 1% 
from the extrapolated values for S > 0. Note that when 
we increase the crystal field S these values change contin- 
uously up to x 9 (l,0) - 0.10 and x 9 (2,0) - 0.28 at the 
multi-critical point. This scenario is quite different from 
that of the Blume-Capel model, which is not a surprise 
since both models are in different universality classes of 
critical behavior. However it is also distinct from the 
scenario of the dilute 4-states Potts model, where the 
scaling dimensions along the critical line are believed to 
be the same as that of the pure 4-states Potts model. If 
the dilution had the same role in both models, we should 
expect for the dilute 4-states Potts a continuous line of 
fixed point, however this was not founds 

IV. SUMMARY AND CONCLUSION 

In this paper we have calculated the phase diagram and 
critical properties of the spin-1 Baxter- Wu model in the 
presence of a crystal field. Our results, based on renor- 
malization group, finite-size scaling and conformal invari- 



ance, show a second-order transition line separated from 
a first-order transition line by a multi-critical point. This 
scenario is in disagreement with that of a previous paper 
by Kinzel et al.jVL where the second-order transition line 
appears only in the limiting case A — > — oo. The crit- 
ical behavior was determined by renormalization group 
and conformal invariance. Despite the phenomenological 
renormalization group be not so conclusive regarding the 
critical exponents, the conformal invariance results indi- 
cate that along the critical line, and even at the multi- 
critical point, the conformal anomaly is the same as that 
of the Baxter-Wu model or the 4-states Potts model, i.e., 
c = 1, in agreement with the scenario expected for the 
dilute 4-states Potts model. However, our results indi- 
cate that the scaling dimensions vary continuously with 
the crystal field. This is an unexpected behavior since 
the reported results for the dilute 4-states Potts model 
indicate a constancy of the scaling dimensions along the 
phase transition line^ and it is expected that both mod- 
els belong to the same universality class of critical be- 
havior. This result implies that either contrary to what 
occurs with thermal and magnetic perturbations the ef- 
fect of dilution is distinct in the Baxter-Wu and in the 
4-states Potts models, or the scenario based in renormal- 
ization group for the dilute 4-states Potts model is wrong. 
It would be interesting to verify our results using differ- 
ent numerical techniques, like Monte Carlo methods. In 
fact, a Monte Carlo simulation for the case 5 — has al- 
ready been done and different exponents as those of the 
pure Baxter-Wu model have been achieved^ Extensive 
Monte Carlo simulations for S ^ Q are in progress and 
will be present elsewhere^ 
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